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A novel lasing scheme for terahertz quantum cascade lasers, based on consecutive phonon-photon-phonon 
emissions per module, is proposed and experimentally demonstrated. The charge transport of the proposed 
structure is modeled using a rate equation formalism. An optimization code based on a genetic algorithm was 
developed to find a four-well design in the GaAs/Alo.25Gao.75As material system that maximizes the product 
of population inversion and oscillator strength at 150 K. The fabricated devices using Au double-metal 
waveguides show lasing at 3.2 THz up to 138 K. The electrical characteristics display no sign of differential 
resistance drop at lasing threshold, which suggests - thanks to the rate equation model - a slow depopulation 
rate of the lower lasing state, a hypothesis confirmed by non-equilibrium Green's function calculations. 
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I. INTRODUCTION 

Nearly a decade after the first demonstration of 
terahertz (THz) quantum cascade lasers (QCL)i, the 
maximum operating temperature (T ma x) of these de- 
vices has reached 199 K, using a three-well resonant 
phonon design 2 -. The high temperature THz-QCLs 
(Tmax <^ 175 K) are mostly designed using resonant 
tunneling based injection and extraction of carriers from 
the lasing states in a GaAs/Alo.isGao.ssAs material 
system 2 - - — . Despite the relatively slow progress on im- 
proving the T max , research efforts remain very intense to 
bring THz QCL into a temperature range achievable with 
thermo-electric coolers. This would open up a profusion 
of THz applications in many areas, including high speed 
communications, pharmacology, non-invasive cross sec- 
tional imaging, quality control, gas and pollution sensing, 
biochemical label-free sensing and security screening^. 

Several theoretical models have been employed to un- 
derstand details of charge transport and optical gain in 
THz QCLs, including density matrix formalism- ~ 9 , non- 
equilibrium Green's function (NEGF)i2r— , and Monte 
Carlo technique a 13 ' 14 . These models have identified sev- 
eral limitations of THz QCLs based on resonant tunnel- 
ing injection (RT-QCL), which make it difficult for the 
lasers with photon energy hu to operate at temperatures 
T > Huj/kB (ks being Boltzmann constant). It is there- 
fore necessary to find alternative designs with novel lasing 
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schemes and/or investigate materials with lower electron 
effective mass^. In this paper we focus on the first ap- 
proach. 

YasudaiS, Kubisil, Kumar— have also pointed out 
the limitations of RT-QCL and proposed a scheme for- 
merly demonstrated in mid-infrared devices^, for THz 
QCLs. In this scheme, carrier injection to the upper las- 
ing state (ULS, also called level 2) is assisted by longi- 
tudinal optical (LO) phonon scattering, and is, hence, 
so-called scattering-assisted QCL (SA-QCL). This type 
of lasing scheme has recently proven its capability to 
significantly surpass the empirical T max ~ hu/kB lim- 
itation of RT-QCL: Kumar et al. have demonstrated a 
four-well GaAs/Al .i5Ga .85As SA-QCL at 1.8 THz with 
T max = 163 K ~ 1.9hw/kB^- A similar design of SA- 
QCL consisting of five wells per module was also demon- 
strated at 4 THz and below 50 K in the lattice-matched 
InGaAs-AlInAs-InP material system 20 . In the following, 
we propose a novel four-well SA-QCL, where the lower 
lasing state (LLS, also called level 1 ) is depopulated with- 
out the mediation of resonant tunneling. The proposed 
design involves only four fairly confined states that are 
isolated from higher energy levels. 

We recall the main drawbacks faced by the RT-QCL. 
The population inversion (AN / N tot ) of three- well reso- 
nant phonon based active regions is limited to 50%. This 
is due to the existence of the injector level, i, where carri- 
ers wait to get resonantly injected into the ULS. Ideally, 
in the coherent transport regime, the ULS holds as many 
carriers as the injector level, which is only half of the 
total available carriers. In a realistic case, with rather 
thick injection barrier and presence of scattering chan- 
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FIG. 1. (Color online) Schematic diagram of a scattering as- 
sisted QCL active region based on a "phonon-photon-phonon" 
configuration. 



nels, the population inversion falls below 50%. More- 
over, as the injector state is populated, thermal backfill- 
ing to the lower lasing state (LLS) cannot be ignored 
at high temperature o 16 ' 21 . It was also predicted that 
the presence of injection and extraction tunneling cou- 
plings broadens the gain with a complicated electric field 
dependency 8 ^ 22 ' 23 , which makes it hard to predict the 
laser frequency. In order to inject carriers selectively to 
the ULS and not the LLS, the injection barrier should 
be chosen meticulously: thick enough to avoid the coun- 
terproductive (also called "wrong") injection of carriers 
to the LLS as well as to reduce the negative differential 
resistance (NDR) at the i-LLS level alignment, and at 
the same time, thin enough to increase the laser's cur- 
rent density dynamic range. It is important to note that 
due to the small ULS and LLS energy spacing (small 
photon energy), the j-ULS and i-LLS level alignments 
happen at similar electric fields. The diagonality of the 
lasing transition was used as a design tool to decouple 
the intermediate tunneling resonances and the injection 
resonance, which occur before and at the design elec- 
tric field respectiviely^i. The tradeoff in designing the 
injection tunneling barrier is even harder to satisfy for 
devices with smaller photon energy. As a consequence, 
the transport across the injection barrier may not be- 
come completely coherent, particularly when considering 
all existing scattering processes, and therefore the peak 
current becomes sensitive to variations of dephasing scat- 
tering during tunneling 18 . 

Refs. [IB, H3 and HH provide comprehensive discussions 
of how SA-QCL designs address the shortcomings of the 
RT-QCL. Figure [Tj illustrates schematically the operation 
principle of the SA-QCLs, with consecutive "phonon- 
photon-phonon" emissions per module. In energy as- 
cending order, the states are named e (for extraction 



from LLS), 1, 2 and i (for injection to ULS). Such de- 
sign schemes rely on an injection energy state, i, lying 
about one LO-phonon energy above the ULS. The car- 
riers are injected from the upstream extraction level, e', 
to the next injection level i by resonant tunneling. Since 
the electric field at e'-i alignment is much larger than 
that for e'-2 and e'-l alignments, the injection reso- 
nance e'-i is well decoupled from the last two resonances 
and a large coupling strength, f2 e ,-, can be used in the 
design. Therefore, it is possible to have thin tunneling 
barriers, while preserving the pumping selectivity to the 
ULS, provided the diagonality between the lasing states 
is adjusted to minimize the relaxation from i to 1. As 
a result, this type of design relies on a large diagonal- 
ity between lasing states, hence the occurrence of pho- 
ton emission assisted tunneling through a rather thick 
"radiative barrier"—. By choosing a thin injection bar- 
rier, one can ensure coherent transport and subsequent 
accumulation of carriers on the state with longest life- 
time - preferably the ULS. This would also leave a small 
concentration of carriers on the short lifetime injection 
state i and extractor state e ', which reduces the effect of 
backfilling to state 1 '. The coherent transport through 
the barrier, additionally, makes current and population 
inversion relatively insensitive to variations of dephas- 
ing processes. As shown in Figure [U for depopulating 
the LLS we simply propose to rely on LO-phonon emis- 
sion scattering, an extraction scheme similar to that of 
Refs. and[25| in two- well RT-QCLs. Our design is sim- 
ilar to the one proposed in Ref. LLji Alternatively, one 
could use the mediation of tunneling to an excited state in 
resonance with LLS before LO-phonon emission scatter- 
ing to the extraction state, a configuration that could be 
named "phonon-photon-tunnel-phonon" 18 i 20 . One could 
also skip the depopulation by LO-phonon emission and 
make the LLS acting as a state that feeds the injection 
state of the next module by resonant tunneling, an ex- 
traction controlled "phonon-photon" configuration which 
offers the advantage of a lower voltage defect^. 

In this work, we focus on "phonon-photon-phonon" 
scheme and evaluate the lasing states (and gain profile) 
that are not perturbed by injection and/or extraction 
tunnelings. In such a scheme, the relative diagonality 
between the four states should be carefully adjusted to 
insure fast injection to ULS and extraction from LLS, 
while minimizing the wrong injection i— >_Z and extrac- 
tion £— s-e channels and at the same time maintaining 
a sufficient oscillator strength between the lasing states. 
In the next section, we will discuss the key figures of 
merit, derived from rate equations and used for design 
optimization. 



II. RATE EQUATIONS 

The key figures of merit, such as population inversion 
and transit time of a "phonon-photon-phonon" design, 
can be derived analytically from rate equations, thanks 



3 



to the operational simplicity of such SA-QCL designs. 
We denote the scattering rate from state I to m as 
gim and the inverse of gi rn as Ti m . The tunneling 
time across the injection barrier can be written as 
Ttun = (1 + A 2 r|)/2f2 2 T||, where HA and Ml are, 
respectively, the detuning energy and tunneling coupling 
strength between the states e and i, and rii denotes 
the phase coherence time constant 2 ^. For the rate 
equation model, the basis of states is taken from the 
eigenstates of "isolated modules", as if the injection 
barrier were very thick. Simplified analytical forms can 
be obtained when backscattering is neglected (i.e. low 
temperature): g m i <C gi m , where E\ > E m . The total 
down scattering rates from levels i and 2 are defined as 
9i = 9i2 + gn + 9ie and g 2 = 321 + ff2e, respectively. For 
the sake of simplicity, we define g\ — g\ e . In the absence 
of stimulated emission, i.e. below threshold (J < J t h), 
the population inversion normalized to the total number 
of carriers per period reads 



Ap = 



T2eS 



n.u 



(1) 



where r 2c ff is called the modified effective lifetime of 
ULS and is defined as 



72 off 



3*2 (9l ~ 521 ) - 32^1 



9i92gi 



(2) 



and r t j is the transit time below threshold defined as 



and Ap w (r 2 i— r-i) / (T tU n+ r u )j respectively. Iftheinter- 
subband relaxation time between the lasing states, t 21 , is 
much longer than the injection and extraction lifetimes, 
and if the tunneling time is much shorter than the tran- 
sit time, the population inversion could approach 100%. 
Eq. [5] also suggests that to avoid accumulation of carri- 
ers behind the tunneling barrier, the transport should be 
coherent (r tun -C t% or equivalently, 2Q, 2 T\\Ti 1). 

In the presence of stimulated emission (J > Jth), the 
population inversion is clamped at Ap tn , which is fixed 
by cavity loss and the gain cross-section. One can show 
that the current density above threshold is 



qN t 



1 + AR th /R 



ih 



Ttu 



(7) 



where t^. is the transit time above threshold and is ap- 
proximately 4ti, assuming efficient and similar injection 
and depopulation rates [ga ~ <?i). More precisely, the 
transit time above threshold becomes 



2 n 1 



9i2 + 9n 

91 + 92e 



(8) 



The fractional change in differential resistance at thresh- 
old, ARth/Rth, is also calculated as 



AR 



th 



R 



th 



R th 



Ap. 



'th" 



9l ~ 92e 
'gi +92e 



(9) 



5*2(31 +321) + .92(231 +9n) 



9i929i 



(3) 



The normalized populations on states i and e can be 
written simply as 



pa 



and 



7"tun ~r 7"tr 



7"tun ~l~ 7~tr 



(4) 



(5) 



A discontinuity of differential resistance is expected 
at threshold, as the current below threshold depends 
on the intersubband relaxation between lasing states, 
T21 (see Eq. [2]), and above threshold this relaxation is 
governed by stimulated emission. From Eqs. [TJ [5] and 
the current density dynamic range gives a condensed 
form of 

Jmax _ 1 + Ai?th/-Rth 2f2 2 T||T 2off 



lim 

coh.trans. \ J^h 



Apth 1 
1 + Ai? th /i? tn f 2cff 



20V||T> 



Ap, 



th 



(10) 



respectively. The current density below threshold 
was estimated in the middle of the tunneling barrier. 
Assuming this layer to be crossed only by tunneling 
currents, the current density can be written simply as 



T< M Pee - Pvi N, 

= qN s = q 



(6) 



where N s is the two-dimensional carrier density per mod- 
ule. In the ideal case, where leakages gn, gi e and 3 2e are 
negligible, the transit time and population inversion con- 
verge towards intuitive expressions: « 2r^ 2 + T21 + T\ 



The second expression shows the importance of using co- 
herent transport across the tunneling barrier in order to 
maximize the dynamic range, a condition reached when 
2f2 2 T||T t ^ ^> 1, or approximatively when 8fl 2 7-|| J> 1. 
The differential internal efficiency of the laser is obtained 



as 



3*2(31 - 321) ~ 323»i 

9i(gi+g2e) 



(11) 



If we define the effective lifetime of level 2 by the stan- 
dard expression, T 2c ff = t 2 (1 — ti/t 2 i), more intuitive 
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expressions for the modified effective lifetime and the 
internal efficiency can be obtained as, 



~ 9i2 9il 

T2e« = T2off T\ and 

9% 9i 



V = 



T2cft 



t-2 off + n 



(12) 
(13) 



respectively. Eg. [T2l shows how the effective ULS lifetime 
is modified by the leakages channels and i— >e that 
are shunting the "correct" trajectory of carriers inside 
a module from i— >2— >1 — >e. The fractional change in 
differential resistance at threshold is closely related to 
the internal efficiency; it can be re- written versus T2 e ff 
as 



th 



AR, 

Rth 



T~2 off + Tl ' 



/0th" 



7"2eff + Tl 



(14) 
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T M„ t (L°) = 0.25 
E = 34.2 meV 
E = 27.7 meV 
S2= 1.14 meV 



21 kV/cm 



an expression which shows similarities with Eq. 1 1 31 From 
Eq. |9l it appears that the discontinuity in differential 
resistance would vanish if the depopulation rate, g\, is 
inefficient and close to the wrong depopulation channel, 
.9i ~ <?2e- Nevertheless, in such adverse conditions, popu- 
lation inversion is still possible, if the transition rate be- 
tween lasing state, 321, is lower than 321 lim = ff2e ~pf^ ■ 
In this situation, the modified effective lifetime, T2 C s, 
would be small ~ (52mm - 52i)/2g 2 ! e (f° r 921 < 521 lim), 
implying the efficiency of the laser will be impeded. Of 
course, the transit time above threshold would be much 
longer than 4rj. In the next section, we describe how the 
rate equation model was employed to design the first iter- 
ation of SA-QCL based on the "phonon-photon-phonon" 
scheme. 



III. DESIGN OF STRUCTURE 

For the first iteration of the "phonon-photon-phonon" 
design, we optimized a figure of merit defined as the 
product of population inversion, oscillator strength and 
inverse of the superperiod length. We chose to maxi- 
mize this figure of merit using four quantum wells in the 
GaAs/Alo.25Gao.75As material system at 21 kV/cm and 
a lattice temperature of 150 K. We assumed a Maxwell- 
Boltzmann distribution for the carriers in all subbands, 
with a characteristic temperature of T e — 200 K. For 
the first iteration, we assumed that transport is domi- 
nated by LO-phonon scattering, and neglected electron- 
impurity and interface roughness (IR) scattering, which 
can be a serious shortcoming. The calculation included 
both forward and backward scattering channels. A more 
accurate figure of merit would include several scattering 
potentials acting on th population rate and dephasing of 
the lasing transition. To avoid leakage channels to higher 
energy states, we deliberately opted for a rather thick in- 
jection barrier of 44 A. Within this model, a genetic 



FIG. 2. (Color online) Conduction band diagram and the 
moduli squared of wavefunctions at 21 kV/cm of the designed 
THz QCL active region based on "phonon-photon-phonon" 
scheme. The "+" signs denote the ionized impurities used for 
doping. The LO-phonon emission assisted scattering times 
r im°' are gi ven f° r the initial kinetic energy Elo — film- 



algorithm based optimization approach was used to find 
the optimum thicknesses of four quantum wells and three 
barriers (the injection tunneling barrier thickness being 
fixed). The energy spacing between the four levels were 
unbound during the optimization process. 

The optimization resulted in the solution schemat- 
ically shown in Figure at 21 kV/cm. Start- 
ing from the injection barrier, it consists of four 
wells and barriers with the layer thicknesses of 
44/62.5/10.9/66.5/22.8/84.8/9.1/61 A - the barriers 
are indicated in bold fonts. The injection barrier is delta- 
doped with Si to 3.25 x 10 10 cm~ 2 , at the center. Inter- 
estingly, with the discussed figure of merit, which only 
includes LO-phonon scattering and the assumption of 
the same electronic temperature of 200 K in all sub- 
bands, the genetic algorithm converged towards a low 
energy spacing between the LLS and the extraction level, 
Ei—E e = 27.7 meV, i.e. 9 meV lower than th bulk GaAs 
LO-phonon energy. This suggests that, within our simpli- 
fied model, the depopulation rate is thermally activated 
(with a 9 meV activation energy), which impedes pop- 
ulation inversion at low temperature. The electrostatic 
energy per module is 76 meV that is approximately twice 
the phonon energy, as proposed by Kubis et ah, to favor 
carrier thermalization within one module^. With the 
choice of thick injection barrier, the tunneling coupling 
strength is rather weak (fif2=1.14 meV), comparing with 
1.5 meV from the first demonstrated THz SA-QCL in 
Ref . OJ. 

In the conduction band diagram shown in Figure [5] at 
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TABLE I. Scattering lifetimes (including LO-phonon, IR and 
e-impurity) in ps at 20 K and 150 K lattice temperatures. 
The IR scattering times are calculated with a correlation 
length of A.ir — 6.5 nm and a mean roughness height of 
Air = 2.83 A. In the fourth column, e-impurity scattering 
times are reported. In the last column the total backscatter- 
ing lifetimes at 150 K are reported. Values are not reported 
for very long n m > 400 ps. In all subbands, electrons were 
distributed according to Maxwell-Boltzmann statistics with a 
temperature, T e , 50 K higher than lattice. 
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21 kV/cm, the four main wavefunctions are completely 
isolated from higher energy states. For simplicity, when 
the wavefunctions are solved for several modules (oppo- 
site case to "isolated modules") we keep the already- 
discussed notation of states e, 1, 2 and i, even when 
the wavefunctions are hybridized from a pair of states 
that are in resonance, like levels e and i" in Figure [2] 
The optimization converged towards excellent wavefunc- 
tion overlaps between i-2 and 1-e, which led to short LO- 
phonon emission scattering times of 0.25 ps for both tran- 
sitions at initial kinetic energy ELo — -E^.ie- The thick 
"radiative barrier" (22.8 A) and the strong diagonality 
between the lasing states lead to an oscillator strength of 
0.39 at 3.46 THz and a LO-phonon emission scattering 
time t^°^ ~ 1.1 ps at 22.4 meV initial kinetic energy 
(-Elo ~ E21 )■ Since the four wavefunctions are fairly di- 
agonal to each other and isolated from higher subbands, 
the free carrier absorption of the active region is probably 
small 28 . This type of design would need to be experimen- 
tally optimized versus the doping density. Without the 
concern of free carrier absorption, the doping could be 
increased; for instance electron-electron scattering could 
favor carrier thermalization and potentially improve de- 
vice performance but on the other hand, e-impurity scat- 
tering could increase the gain linewidth, particularly at 
high temperatures^. Even though the optimization pro- 
cess was performed with only LO-phonon scattering, the 
subsequent simulations include LO-phonon, IR and e- 
impurity scattering potentials. For the sake of compar- 
ison, the scattering lifetimes of the six channels, shown 
in Figure [1] are listed in Table fl] As expected, with the 
assumption of Maxwell-Boltzmann carrier distributions, 
the total scattering time for depopulation, T le , is shorter 
at elevated temperatures, which is due to the 9 meV ac- 
tivation energy of LO-phonon scattering. The IR scat- 
tering is expected to give a large contribution between 
vertical states (i.e. i-2 and i-e), but we find it is also not 
negligible between the two diagonal lasing states^. The 




Electric Field (kV/cm) 



FIG. 3. (Color online) Results of simulations based on a rate 
equation model performed on the structure proposed in Fig- 
ure [2] a) Different characteristic times at 20 K (thick lines) 
and 150 K (thin lines)- 70 K and 200 K for electrons, respec- 
tively. As discussed in section [II] T tun is tunneling time (solid 
line); t£ (dot line) and t£ (dash dot line) are transit times 
before and after threshold; r,: is injection state lifetime (dash 
line); T2 E a is standard effective lifetime of ULS (short dash 
line); and T2eS is modified effective lifetime (dash dot line), 
b) Normalized populations of the four states at 20 and 150 
K lattice temperatures, c) Current density, lasing frequency 
and optical gain versus electric field at 20 and 150 K lattice 
temperatures. 



IR scattering contributes 76% (44%) to the total scatter- 
ing rate 1721 at 20 K (150 K). In the next optimization 
iteration, it would be important to include IR scattering 
in the figure of merit. Several wafers based on this simple 
design with different growth conditions could be tested 
in order to elucidate the role of IR scattering in THz 
QCL, a subject which is debated in the litteratur o 11 ' 12 ' 29 . 
A more elaborate figure of merit could also include the 
gain bandwidth, estimated by various scattering poten- 
tials (IR scattering is likely to contribute the mos t 30 ' 31 ). 
These refined models would adjust the relative diagonal- 
ity between the four states. Table U shows that electron- 
impurity scattering gives a small contribution. 

The current density, optical gain, populations and dif- 
ferent characteristic lifetimes were calculated within the 
rate equation model. The results versus electric field at 
two temperatures (20 K and 150 K) are shown in Fig- 
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FIG. 4. (Color online) Left axis: peak gain at 21 kV/cm 
versus the electronic temperature from simulations based on 
rate equation model. Right axis: intersubband recombination 
time between lasing states, T21, and depopulation time con- 
stant, rie, the difference T21 — ru, and the modified effective 
lifetime of ULS, en- 



ure [3] In the simulation, the coherence time constant, 
T||, has a pure dephasing time component, r* = 0.35 ps, 
which is assumed to be temperature independent for sim- 
plicity. The optical gain was estimated by assuming a 
constant intersubband linewidth of 1 THz. The different 
lifetimes, as introduced in section UH are plotted in Fig- 
ure [3]-a. We can see the tunneling time, r tun is slightly 
shorter than the transit time above threshold, r t ^ (by a 
factor » 2.2 at 20 K), meaning that the current dynamic 
range is not optimized (see Eq. 1101) . The tunneling time 
is substantially longer than the injection state lifetime, 
T,-, by a factor of 2.5, which indicates transport is not 
coherent through the barrier. As a result, at the design 
electric field, the populations on the states e and i are 
very different. In other words, carriers tend to accumu- 
late on the extraction state, and hence promote backfill- 
ing to the LLS (see Eq. [5]) . Figure [3}b illustrates the 
effect of backfilling on populations pn, p ee and P22, par- 
ticularly at electric fields below 21 kV/cm. At the design 
electric field, the effect of backfilling at high temperature 
is mitigated by the faster relaxation between laser states 
T21 and depopulation T± e . At the design electric field, the 
population of level i increases with temperature, as the 
transit time t£ for carriers contributing to current de- 
creases, while its lifetime remains unchanged. As shown 
in Figure [3}c, the slight decrease in current density at 
higher temperature is related to backfilling to the LLS, 
meaning that less carriers participate in the transport. 

The peak gain value at 21 kV/cm was simulated versus 
temperature and the results are displayed in Figure 3) 
According to the model presented in section [TlJ a gain 
peak value of more than 40 cm -1 is maintained up to 
T e — 230 K. The surprising rise of gain from low temper- 
ature up to T e — 110 K originates from the assumption 
of Maxwell-Boltzmann distributions. This is due to the 
temperature-induced enhancement of the depopulation 



rate (1— >e) being greater than that of the intersubband 
relaxation between the lasing states (2— ll), at low tem- 
peratures. Indeed, the latter has an activation energy 
of 22.4 meV by LO-phonon scattering (versus 9 meV for 
<?i e ) and involves a pair of states that are more diago- 
nal than the pair 1-e. Consequently at low temperature, 
321 has a stronger contribution from IR scattering than 
gi e , which is harder to change by LO-phonon scattering 
(see Table U and right axis of Figure |4|. We can, further, 
see from the Figure 0] that the temperature dependance 
of T21 — ri e and the modified effective lifetime, T2eSi are 
mirrored in the gain. In the next section we will describe 
the experimental results obtained with the structure pro- 
posed in Figure[5J 



IV. EXPERIMENTAL RESULTS 

The design presented in the Figure [5] was grown on a 
semi-insulating GaAs substrate by molecular beam epi- 
taxy with 276 repeats to obtain a 10 /im thick active re- 
gion. The active region was sandwiched between a bot- 
tom stack 100 nm of 3 x 10 18 cm -3 n + GaAs followed 
by 20 nm of non-intentionally doped GaAs spacer and a 
top stack of 8 x 10 17 /5 x 10 18 cm~ 3 (40 nm/50 nm) n + 
GaAs layers followed by 10 nm of 5 x 10 19 cm~ 3 n + GaAs 
layer and capped with 3 nm of low-temperature grown 
GaAs. The lower doping concentration of the bottom 
contact layer (3 x 10 18 cm -3 ) and the 20 nm spacer were 
meant to limit diffusion of Si dopants to the first period, 
which could induce intermixing between the first layers 
of the active region. Indeed, by transmission electron mi- 
croscope imaging, we observed blurred interfaces on the 
first three barriers of the structures that were grown with 
100 nm of 5 x 10 18 cm~ 3 bottom n + GaAs followed by 
10 nm of GaAs spacer—. Special emphasis was put on 
minimizing the drift of fluxes for Ga and Al during this 
long growth process. The X-ray diffraction rocking curve 
could be fitted perfectly with nominal parameters, with 
no extra broadening of satellites peaks, confirming the 
excellent stability of the growth rates (better than 0.5%), 
throughout the active region. The wafers were processed 
into THz QCL structures with Au double metal waveg- 
uides. Devices have the following dimensions: ~ 144 /im 
wide ridges with ^120 /im wide top Ti/Au metallization 
forming a Schottky contact, and ~ 1 mm long Fabry- 
Perot resonator. An In-Au wafer bonding technique was 
use d 33 ' 34 and the ridges were fabricated by reactive-ion 
etching. The laser bars were indium soldered (epi-layer 
side up) on silicon carriers and then mounted in a He 
closed-cycle cryostat for measurements. 

Figure [5] shows the light(L)-current density (J)-voltage 
(V) characteristics of the fabricated device from 10 K to 
138 K in pulsed mode, with pulse duration of 250 ns and 
repetition rate of 1 kHz. The device showed a T max of 
138 K at a current density of 1.65 kA/cm 2 . At 10 K, the 
threshold current density is 1.17 kA/cm 2 and the maxi- 
mum current density is 1.55 kA/cm 2 resulting in a small 
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FIG. 5. (Color online) Right axis: collected THz light (opti- 
cal output power) versus current density curves for Au double- 
metal THz QCLs based on the "phonon-photon-phonon" las- 
ing scheme at different heat sink temperatures. Left axis: 
voltage (V) versus current density (J) at 10 K and 130 K. The 
dashed horizontal lines highlight the position of resonances in 
the V-J plot. The devices are ~ 144 /xm wide, 1 mm long and 
are fabricated using TiAu metal contacts. The bias is applied 
in pulsed mode (pulse width = 250 ns, repetition rate = 1 
kHz). The top inset shows a zoom of THz light and differ- 
ential resistance versus current density around threshold at 
10 K. The bottom inset shows the threshold density versus 
temperature. 



ratio of J m ax/>/th = 1-325. The laser action starts at 
20.25 V and stops at 21.8 V where a NDR is observed. 
By subtracting the 0.8 V Schottky voltage drop on the 
top contact* 3 ^, the NDR voltage corresponds to the design 
electric field of 21 kV/cm. A zoom- view of the differential 
resistance, R ac , versus current density around threshold 
at 10 K (inset of Figured]) did not show a clear disconti- 
nuity at threshold. Several measurements around thresh- 
old with different voltage steps were attempted but no 
clear change in the slope of the J-V characteristics could 
be detected. Analyzing this observation with the model 
discussed in section|TT]suggests that the depopulation rate 
1— is not as efficient as expected. Moreover, according 
to Eq. [TUJ the ratio J max / Jth can be approximated by 
T^cff/ApthT't'r ■ The vanishing discontinuity of differential 
resistance is consistent with the low output power of the 
laser. It suggests that the effective lifetime of the ULS 
is much shorter than expected, and consequently reduces 
the internal efficiency of the laser (see Eq. IT3| . The slow 
rate of depopulation could be related to the E\ e energy 
spacing, being 9 meV smaller than the GaAs LO-phonon 
energy. As the LO-phonon emission scattering shows 
a strong resonance at initial kinetic energy -Elo ~ E\ ei 
any deviation from the Maxwell-Boltzmann distribution 
could seriously alter the average depopulation rate g\ e - 
a crucial parameter for laser action. If thcrmalization of 
carriers by e-e intrasubband scattering is not significantly 
more efficient than intersubband scattering processes, the 
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FIG. 6. (Color online) Conduction band diagram and the 
moduli squared of wavefunctions at a) 4.4 kV/cm and b) 8.7 
kV/cm. 



carriers accumulate in the subband 1 up to the kinetic 
energy of 9 meV, in contrast to the Fermi-Dirac or the 
Maxwell-Boltzmann distribution o 35 ' 36 . The relative sim- 
plicity of this device would make it a good candidate for 
studies of electron temperature by microprobe photolu- 
minescence experiments 3 ^. 

At 10 K, the voltage versus current density plot shows 
a tunneling resonance with a shoulder at ~5.6 V and a 
strong NDR at 9.4 V (see left axis on Figure [5]). Tak- 
ing into account the 0.8 V Schottky voltage drop on 
the top contact, the electric field of these features cor- 
respond well to the alignment of the extraction state 
e with levels 1 and 2, predicted at 4.4 kV/cm and 8.7 
kV/cm, respectively. Figure |B] depicts the wavefunction 
configurations for these electric fields. The anticross- 
ing of the lasing states occurs at 13 kV/cm, implying 
that below 13 kV/cm level 2(1) is located on the down- 
stream(upstream) side of the module. At 4.4 kV/cm, 
levels e n _i and l n are aligned with a coupling strength 
of fi e i = 0.385 meV, while the next set of anticrossed 
states, e n and l a +i, are 16 meV below (qL£, £ being 
the electric field). It would suggest that the transition 
U\ £n is mainly assisted by IR scattering, with a few ps 
lifetime. At 20 K and for the e — > 1 tunneling transport 
4f2g 1 T|| el T 1 ss 1, suggesting that transport is "in between" 
the coherent and incoherent regimes. Interestingly, at the 
same electric field, levels 2 n +i and in+2 anticross and are 
also aligned with l n and e n _i. These four states could 
form a weakly coupled miniband (see Figure [6]-a). One 
could imagine transport occurring two periods at a time 
from e n -i to 2 n+ i, since 2qL£ = 32 meV w -Elo- There- 
fore, one can propose that the shoulder in V-J plot at 
~5.6 V, which is also slightly higher than e-1 alignment 
voltage (4.4+0.8=5.2 V), corresponds to an electric field 
for which 2qL£ s» Elo (4-8 kV/cm when taking into ac- 
count the Schottky voltage drop). At 8.7 kV/cm, levels 
e n _i and 2 n are aligned with a small coupling strength 
^e2 = 0.24 meV. Transition to the next lower extraction 
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FIG. 7. (Color online) THz spectra recorded for different 
biases and temperatures. The bias per module is reported 
between brackets. 
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FIG. 8. Energy resolved electron density and the most sig- 
nificant basis states |ip a (2;)| 2 (a. u.) at a bias of 78 mV per 
period and a temperature of 140 K. 



state e n should be efficient, as qL£ = 32 meV « -Elo- 
The figure of merit of e — > 2 tunneling transport is 
4f2g 2 T|| e2 T2 « 0.04 < 1 at 20 K, suggesting that this 
transport channel is incoherent, and hence very depen- 
dent on phase coherence time constant&i 8 .. 

Figure [7] depicts the measured spectra for the las- 
ing device at different biases and temperatures. At 
10 K and slightly above threshold, the emission line is 
at ^2.83 THz; at higher bias, it eventually shifts to 
higher frequency, 3.23 THz. According to the solution 
of Schrodinger equation, in the range of electric fields 
for which the spectra were recorded (~19.7-21 kV/cm), 
the frequency spacing between the lasing states should 
move from 3.05 to 3.46 THz (see Figure [3]-c) . At this 
stage, it is difficult to strongly affirm that this differ- 
ence of ~ 0.9 meV in photon energy is an unequivocal 
sign of a physical phenomenon we did take into account. 
However, since we have identified a signature of a sig- 
nificantly long LLS lifetime (no detectable discontinuity 
of i? ac at threshold) and therefore, a residual popula- 
tion on this state, we cannot exclude the occurrence of 
Bloch gain between the lasing state d 38 ' 39 . The disper- 
sive contribution of Bloch gain would drag the peak gain 
frequency to a lower value than z^i- This design would 
be an interesting candidate for gain measurements by 
THz time-domain spectroscop y 40 ' 41 to see if any spectral 
asymmetry in gain can be detected. In the next section, 
the results of NEGF simulations for this QCL structure 
are presented. 



V. SIMULATIONS BY NONEQU I LIBRIUM GREEN'S 
FUNCTION THEORY 

We have simulated current density, gain spectra and 
carrier densities with the method of NEGE 4 ^ -. The 
lesser Green's function G < is related to the density ma- 



trix p a p(k) = —i J H^G^g(k, E). The calculation is done 

with the basis states A^ 1 ^ 2 ^2 a ip a (z)e~ l]lr , where the 
bold typeface signifies two dimensional in-plane vectors 
and A is the cross-section. This gives the energy resolved 
electron density 

n(z,E) = jJ2J2^ G ^( k ' E H( z )^( z )}- ( 15 ) 

k a/3 

The current density (in A/cm 2 ) is 

m = - f d dzj( Z ,t) = 4ee <;w(m),(16) 
aj ° a a p k 

where 

W «^\j/ Z 2^(z) 

( *, s d<p a {z) d^ 3 (z) \ 

[M*)-9- z K-vaWjm 

G < is calculated on the basis of the envelope function 
Hamiltonian of the perfect QCL structure, where pertur- 
bations by interface roughness, ionized impurities, and 
phonons are treated by self-energies in the self-consistent 
Born approximation. The electron-electron interaction is 
taken into account within the mean-field approximation. 
In the simulations we have used an effective roughness 
height of 2 A, and an average roughness distance of 100 
A. The simulated electron densities are shown in Fig- 
ure along with the significant basis states <p a (z). The 
darker regions close to the bottom of the subbands in- 
dicate that carriers are not thcrmalized, contrary to the 
assumption used in the previous simulations based on 
rate equations. We can clearly identify in which states 
and up to which energy carriers accumulate. This graph 
suggests that electrons get trapped in the LLS since the 
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FIG. 9. (Color online) Current density versus voltage bias 
per period: experimental data (scattered points) at 10 K in 
pulsed mode and comparison with simulations by nonequilib- 
rium Green's function formalism. The position of the first 
resonance at ~4.8 kV/cm in the experimental data is high- 
lighted by a vertical arrow and is attributed to sequential 
tunneling from e to 1 and then 2 in three consecutive peri- 
ods. The NEGF simulation are performed for higher temper- 
atures than the experimental data but the simulated results 
are weakly temperature sensitive. Moreover, by the end of 
the current pulses (250 ns) it is well known that the lattice 
temperature can be several tens of kelvin higher than heat 
sink. The large fluctuations around 40-60 meV are artifacts 
of the simulation. The device was driven in current mode 
which explains why the valley of the second NDR predicted 
by theory is not observed. 



energy gap between this level and the injector state is 9 
meV lower than GaAs LO-phonon energy. Moreover, we 
observe a slight accumulation of carriers at the bottom 
of the upper subband since the energy gap between this 
state and the ULS of the next period does not match ex- 
actly the LO-phonon energy. One can see that the shaded 
region associated with the extraction state is darker than 
that of the injection state: in other words, the popula- 
tion of the extraction state is significantly higher than the 
injection state due to the thick injection barrier. We re- 
call the tunneling coupling strength is rather small: 1.14 
meV. This makes the emptying of the extraction level 
slow, hence promoting the backfilling to LLS. 

In Figure [9l we find excellent agreement with the mea- 
sured data for the current densities. The second NDR 
at 30 mV predicted by NEGF simulations is confirmed 
by the experiment. The device was driven in current 
mode hence, on the oscilloscope, the voltage suddenly 
changed from 30 mV to 60 mV. The first tunneling res- 
onance is predicted at around 4.3 kV/cm whereas the 
experimental current density shows a shoulder at 4.8 
kV/cm (see the vertical arrow in Figure which was 
attributed to sequential tunneling across three modules 
from e n _i —> l n — > 2 n+ i in section [IVJ NEGF simu- 
lations performed around the first resonance confirmed 
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FIG. 10. (Color online) Gain spectra simulated using NEGF 
model at T=140 K and different biases per module. 



a non negligible transport channel by phonon emission 
(and more marginally by other scattering sources) from 
the anticrossed levels 2 and i. Experimentally, the cur- 
rent is more phonon driven at the first resonance than 
predicted in NEGF simulations. For this latter model, 
the scattering potentials are very efficient when levels e 
and 1 are best aligned and therefore, give rise to a current 
peaked at 4.3 kV/cm. 

Finally, gain spectra have been simulated using NEGF 
at different temperatures and voltage biases. The results 
of these simulations for a fixed bias of 21.5 kV/cm (78 mV 
per module) and a temperature of 140 K are displayed in 
Figure fTUl It is important to point out that simulations at 
low temperatures showed a maximum gain and peak gain 
frequency lower than expected from the measurements. 
There were numerical problems at low temperature due 
to the lack of e-e scattering in the simulations, whereas 
at 140 K the acoustic phonons assist to a higher degree 
in the thermalization of the carriers. An underestimated 
thermalization of the carriers at low temperature would 
enhance the trapping of electrons on the LLS, and sub- 
sequently result in a reduced population inversion and a 
higher contribution from Bloch gain. At 140 K, a gain 
of 37 cm -1 at 3.4 THz is predicted, a value that is con- 
sistent with previous measurements of waveguide loss by 
cavity frequency pulling^. Theoretically, from 72 to 76 
mV per module, the peak gain frequency moves from 
2.8 to 3.24 THz, while experimentally, the laser emission 
shifts from 2.8 to 3.23 THz in a similar bias range (71.6 
to 76.3 mV) at 10 K. This latter comparison would sup- 
port the hypothesis that the simulated thermalization at 
140 K by NEGF is close to the actual thermalization 
at lower temperatures. At 76 mV per module, the gain 
bandwidth derived by NEGF simulations is ^1.38 THz, 
a value slightly broader than the assumed value of 1 THz, 
used in the simulations based on rate equations. 
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VI. CONCLUSION 



This paper presented a novel THz SA-QCL design with 
consecutive "phonon-photon-phonon" emissions. The 
simplicity of the lasing scheme, made it convenient to 
calculate the energy level populations and lifetimes using 
a rate equation model. The first iteration structure was 
designed by optimizing a figure of merit defined as the 
product of population inversion and oscillator strength. 
The optimum design at 150 K only took LO-phonon scat- 
tering into account in transport calculations and assumed 
a Maxwell-Boltzmann distribution of carriers with a tem- 
perature of 200 K in the subbands. In the designed struc- 
ture, the lasing states were fairly isolated from higher 
bands and not perturbed by tunneling. The fabricated 
device using Au double metal waveguide lased up to 138 
K. The voltage - current density characteristic of the de- 
vice showed no obvious discontinuity in the differential 
resistance at threshold. The rate equation model sug- 
gests that this effect stems from rather slow depopulation 
of the LLS due to smaller-than-phonon-energy E\ e . This 
hypothesis was confirmed by NEGF calculations. Both 
rate equation and NEGF calculations indicated that the 
tunneling barrier is the bottleneck of the carrier trans- 
port and the carriers are piled up behind it. 

Using the models discussed in this paper, we identi- 
fied the shortcomings of the first iteration design and 
addressed them in detail for the next generations of THz 
SA-QCL designs. We found out that it is important to 
maximize the phonon scattering assisted extraction and 
injection rates, by a better optimization of energy spacing 
between corresponding subbands. Moreover, one has to 
take interface roughness scattering into account in the op- 
timization process. It was found that interface roughness 
scattering, particularly between the lasing states, is com- 
parable to the LO-phonon scattering rates and should 
not be neglected. Its effect on population dynamics and 
linewidth of the lasing transition also has to be consid- 
ered for the next generation of THz SA-QCL designs. 

The proposed lasing scheme in this work benefits from 
a rather simple operation principle. Consequently, such a 
4-level quantum mechanical structure provides a unique 
platform to study experimentally various debated effects 
in THz QCLs, including effect of interface roughness on 
population dynamics and gain linewidth. The authors 
suggest conducting a comprehensive doping study to elu- 
cidate the debated topic of free carrier absorption and 
the role of e-e and e-impurity scatterings in THz QCL. 
Furthermore, one can explore the effect of injection tun- 
neling barrier thickness on the coherence of the transport 
and its effect on the device performance. Further under- 
standing of the carrier dynamics and gain is possible by 
microprobe photoluminescence experiments to study the 
electron distributions in the different subbands and by 
THz time-domain spectroscopy to investigate the possi- 
bility of Bloch gain in such a structure, respectively. 
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